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Abstract. We investigate the behaviour of aft dynamos 
with a dynamic a, whose evolution is governed by the im- 
balance between a driving and a damping term. We focus 
on truncated versions of such dynamo models which are 
often studied in connection with solar and stellar variabil- 
ity. Given the approximate nature of such models, it is 
important to study how robust they are with respect to 
reasonable changes in the formulation of the driving and 
damping terms. For each case, we also study the effects 
of changes of the dynamo number and its sign, the trun- 
cation order and initial conditions. Our results show that 
changes in the formulation of the driving term have im- 

. portant consequences for the dynamical behaviour of such 
systems, with the detailed nature of these effects depend- 
ing crucially on the form of the driving term assumed, the 
value and the sign of the dynamo number and the initial 

• conditions. On the other hand, the change in the damping 
term considered here seems to produce little qualitative 
effect. 

Key words: Sun and stars: magnetic fields - mean field 
dynamos - nonlinear dynamics - fragility - chaos 



1. Introduction 

It is commonly believed that the observed solar and stellar 
variabilities have their origin in the hydromagnetic dy- 
namos associated with turbulent convection zones. Nu- 
merical studies have been made using the full magneto- 
hydrodynamical partial differential equations (PDE), 
which reproduce some features of solar and stellar dy- 
namos (e.g. Gilman 1983). Such models are fairly com- 
plex and do not allow extensive parameter surveys. As a 

Send offprint requests to: Eurico Covas 
* e-mail: E.O.Covas@qmw.ac.uk 
** e-mail: A.S.Tworkowski@qmw.ac.uk 
*** e-mail: Axel.Brandenburg@ncl.ac.uk 
1 e-mail: reza@maths.qmw.ac.uk 



result, a number of alternatives to the direct integration 
of PDE have been pursued. Among these has been the 
employment of the mean field dynamo formalism (Krause 
& Radler 1980) in order to construct various types of dy- 
namos, such as afl dynamo models. Despite the fact that 
such models have been shown to be capable of produc- 
ing a large number of observationally relevant modes of 
behaviour, ranging from stationary to chaotic (c.f. Bran- 
denburg et al. 1989a, b; Tavakol et al. 1995), they nev- 
ertheless involve a number of unknown features such as 
the exact nature of the nonlinearities involved. Further- 
more, in order to clarify the origin of dynamical modes 
of behaviour observed in dynamo models, further simpli- 
fications of these models have been considered, involving 
low dimensional truncations of the governing PDE. Such 
models have also been shown to be capable of produc- 
ing a number of important features of stellar variability 
including periodic, intermittent and chaotic modes of be- 
haviour (Zeldovich et al. 1983; Weiss et al. 1984; Feudel 
et al. 1993). 

Now given that these models are cheaper to integrate 
and more transparent to study, it would be very useful 
if we could employ them as diagnostic tools in order to 
study the effects of introducing different parametrisations 
and nonlinearities involved. The problem, however, is that 
these low dimensional models involve severe approxima- 
tions, and therefore in order to be able to take the results 
produced by them as physically relevant, it is important 
that they remain robust under changes which fall within 
the domain of the approximations assumed. This is par- 
ticularly of importance since on the basis of results from 
dynamical systems theory, structurally stable systems are 
not everywhere dense in the space of dynamical systems 
(Smale 1966), in the sense that small changes in mod- 
els can produce qualitatively important changes in their 
dynamics. In this way the appropriate theoretical frame- 
work for the construction of mathematical models and the 
analysis of observational data may turn out to be that of 
structural fragility (Tavakol & Ellis 1988; Coley & Tavakol 
1992; Tavakol et al. 1995). 
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Here as examples of such changes we shall consider first 
changes in the order of truncation and then changes in the 
details of the physics assumed. Regarding the former, a 
number of attempts have already been made to study the 
effects of increasing the truncation order on the resulting 
dynamics. For example, Schmalz & Stix (1991) (hereafter 
referred to as S&S91) have looked at the detailed dynam- 
ics of the low dimensional truncations of the mean field 
dynamo equations and have studied what happens as the 
order of the truncation is increased, while Tobias et al. 
(1995) have employed normal form theory to construct a 
robust minimal third order model which exhibits both the 
modulation of basic cycles and chaos. These studies have 
shown that low dimensional models can capture a number 
of important dynamical features of the dynamo models. 

Our aim in this paper is complementary to that of the 
above authors. We take a detailed look at the results in 
S&S91 and ask to what extent these results remain ro- 
bust as reasonable changes are made to the details of the 
physics employed, and in each case we study how such 
changes affect the dynamical behaviour of different trun- 
cations. 

2. Models with dynamical a 

The starting point of the truncated dynamical a models 
considered in S&S91 is the mean field induction equation 



dB 

~dt 



= V x (v x B + aB -r] t V x B), 



(1) 



where B and v are the mean magnetic field and the mean 
velocity, respectively. The turbulent magnetic diffusitiv- 
ity rjt and the coefficient a, which relates the mean elec- 
trical current arising in helical turbulence (the a-effect) 
to the mean magnetic field, both arise from the correla- 
tion of small scale (turbulent) velocity and magnetic fields 
(Krause & Radler 1980). 

S&S91 employ an axisymmetrical configuration with 
one spatial dimension x, which corresponds to a latitude 
coordinate and a longitudinal velocity with a constant ra- 
dial gradient (the vertical shear uj )- The magnetic field 
takes the form 



B 



0,B 6 



1 dA 4 
' R dx 



(2) 



where is the ^-component (latitudinal) of the mag- 
netic vector potential, B<j, the ^-component of B and x is 
measured in terms of the stellar radius R. These assump- 
tions allow Eq. (1) to be split into 



dA 4 
dt 

mt 

dt 



Vt d 2 A^ 
R 2 dx 2 
Vt d 2 B <p 



+ 



aB^, 
u dA 4 



R 2 dx 2 R dx 



(3) 
(4) 



In S&S91, a is divided into a static (kinematic) and a 
dynamic (magnetic) part: a — aocosx — Qm(^), with its 



time-dependent part ajw(i) satisfying an evolution equa- 
tion in the form 



daM 
dt 



A(a M ) + f(B), 



(5) 



where A is a damping operator and f(B) is a pseudo- 
scalar that is quadratic in the magnetic filed. 

It has been argued that the a effect is quenched by the 
current helicity density J B, which in turn is governed by 
a dynamical equation (Klecorin & Ruzmaikin 1982; Zel- 
dovich et al. 1983). The reason the feedback (quenching) 
is not instantaneous is a consequence of the fact that the 
magnetic helicity is conserved in the absence of diffusion 
or boundary effects. Such models have been investigated 
recently by Kleeorin et al. (1995). In S&S91 a truncated 
version of yet another model was studied, in which instead 
of the current helicity density, the magnetic helicity den- 
sity A B, or rather A^B^,, was used. Their model was 
motivated on heuristic grounds. Bifurcation properties of 
a truncated version of a similar model, but with a dif- 
ferent damping term, have been studied by Feudel et al. 
(1993). Our present investigation is thus motivated par- 
tially by the variety of models presented in the literature. 
It is important to know what is the effect of the dynamical 
feedback and how the different representations affect the 
results. 

To proceed S&S91 specify the feedback in the following 

way 



f(B) cx A^B^ 



(6) 



and then look at various N— modal truncations of these 
equations and study what happens to the dynamical be- 
haviour of the resulting systems as N is increased. 

To do this it is convenient to transform these equations 
into a non-dimensional form. This can be done by employ- 
ing a reference field Bq, measuring time in units of R 2 /vt 
and defining the following non-dimensional quantities 



A 



Ru 
BoVt 



A, 



Vt 



B=^t. 
Bo 

Vl ' 



c 



(XmB?LOq 



(7) 



where v t is the turbulent diffusivity. Equations (3), (4) 
and (5) with the damping operator taken to be 



A = 



R 2 dx 2 ' 



(8) 



can then be rewritten in the following non-dimensional 
forms: 



d 2 A 



OA _ 

dt dx 2 

dB d 2 B OA 

dC d 2 C 
~dt 



— — = —r + DBcosx-CB, 



dx 2 ^ dx ' 



dx 2 



+ AB. 



(9) 
(10) 
(11) 
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Now considering the interval < x < n (which corre- 
sponds to the full range of latitudes), taking the bound- 
ary conditions at x = and x = tt to be given by 
A = B = C = and using a spectral expansion of the 
form 



iV 



A = A n (t) sin ra, 

n=l 
N 

B = B n (t) sin hi, 

n=l 
N 

C = ^C„(t) sin 



nx, 



(12) 
(13) 
(14) 



n=l 



allows the set of Eqs. (9-11) to be transformed into the 
form 



^ = -n 2 A n + j(B n -! + B n+1 ) 

N N 

+ ^F(n,m,l)B m Ci, 



(15) 



m=l 1=1 



N 



OB 

-n 2 B n + ^ G(n,m)A m , 

m— 1 



-vn 



dt 

dC n 
dt 

where 

F(n, to, I) = 



N N 



: c n - j2J2 F ( n > m > l ) A ™Bi 



m=l 1=1 



8nml 



(16) 



(17) 



(18) 



ir(n + to + l)(n + to — l){n — to + l)(n — to — I) ' 
if n + m + I is odd and F(n 7 m, I) = otherwise and 
4uto 



G(n, to) 



7r(n 2 — to 2 ) ' 



(19) 



if n + m is odd and G(n, to) = otherwise. 

These rules enable the system to describe fields which 
are strictly symmetric (i.e. having only components B n 
with odd n and A n and C n with even n) or strictly an- 
tisymmetric (i.e. having only components A n with odd n 
and B n and C„ with even n) with respect to a; = tt/2, pro- 
vided the initial conditions have either of these parities. 

Using these equations, S&S91 studied a number of such 
truncations numerically by varying the dynamo number D 
at each truncation N. Their main conclusions were: 

1. With the choice of the driving term / given by Eq. 
(6) the antisymmetric truncation with the smallest 
non-trivial indices is identical with the Lorenz system 
(Lorcnz 1963). 

2. Different truncations are capable of producing station- 
ary, oscillatory and chaotic modes of behaviour. They 
also make observations about the changes in the route 



to chaos, and conclude that, as N is increased, the 
route changes from period doubling to the Rucllc 
Takens-Newhouse scenario (Ruelle & Takens 1971; 
Ncwhouse et al. 1978). 

3. The qualitative behaviour of the truncations stabilises 
as the number of modes is increased and in particular 
for TV > 6. As an example they observe that as N 
is increased the limit cycles remain stable for larger 
dynamo numbers. 

4. They also discuss very briefly the D < case, observ- 
ing that the N = 2 case is always a stable fixed point 
and that for N > 6 the antisymmetric limit cycle be- 
comes unstable via a saddle node bifurcation 1 . 

Now, as mentioned above, there are arguments in sup- 
port of both the form of the driving term as well as the 
damping term being different (Kleeorin et al. 1995). So 
as a first step, we shall study, in the next section, how 
robust the results in S&S91 are with respect to various 
physically justified changes in the driving term that have 
been considered in the literature in Eq. (5). In Sect. 4 we 
study the effects of changes in the damping term. 

3. Robustness with respect to changes in the driv- 
ing term 

The general physically motivated choice for the driving 
term is given by Kleeorin & Ruzmaikin (1982), Zeldovich 
et al. (1983) and Kleeorin et al. (1995) to be in the form 



/ = Wi J ■ B + W 2 a\B\ 



(20) 



where W\ and W2 are constants. To study the effects of 
each term separately, we shall proceed by considering the 
cases Wi ^ {W 2 = 0) and W x ^ (W 2 ^ 0) in the 
following sections. 

3.1. Case (I): f = W X J ■ B + W 2 a\B\ 2 , with W x ^ 
(W 2 = 0) 

Taking / to be of the form / oc J ■ B, substituting for B 
from Eq. (2) and recalling that J = V x B we obtain 



J • B = (V x B) ■ B 



B 2 T) t fdAdB d 2 A 



irriTr-io-- 8 ' ( 21 ) 



R 3 uiq \ dx dx dx 2 
which allows Eq. (11) to be written as 



dC _ v d 2 C OAdB d 2 A B 
dt dx 2 dx dx dx 2 



(22) 



Proceeding in a similar way as in previous section we 
obtain an identical set of differential equations to those 



1 Care must be taken when speaking of antisymmetric solu- 
tions. In our studies we mean strictly antisymmetric solutions, 
while in S&S91 these also refers to the antisymmetric part of 
mixed parity solutions. 
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obtained in S&S91, except that Eq. (17) is now changed 
to 



dC n _ 

at 

where 



N N 



m=l 1=1 



H(n, m, I) 
4 



(23) 



(24) 



nml(-n 2 + 3m 2 + I 2 



ir (n + I + m){n + I — m){n — I + m)(n — I — m) 1 

if n + m + I is odd and H(n, m,l) = otherwise. The 
function H is clearly different from F unless F = 0, in 
which case H is also equal to zero. 

For this system we can study also the pure antisym- 
metric and symmetric solutions, but for the sake of com- 
parison with the results in S&S91 we confined ourselves 
to the antisymmetric solutions. 

Now for the case of N = 2, the Eqs. (15), (16) and 
(23) become 



dA 1 

dt 
dB 2 

dt 
dC 2 

dt 



-A 1 + 

-4B 2 - 



DB 2 32B 2 C 2 



= -4z^C 2 + 



2 
16 



15tt 



5ir 



which upon using the transformations 



64 



32 



64 



(25) 
(26) 
(27) 

(28) 



result, as in S&S91, in the usual Lorenz equations (Lorenz 
1963), with the control parameters given by u = 4, b = &v, 
and r = D/3tt. To be compatible with S&S91 we also used 
v = 0.5 throughout 2 . 

Since our aim is to study the qualitative effects brought 
about by the changes in the form of /, we will not 
delve deeply into the details of the dynamics, such as the 
routes to chaos, and concentrate instead on the occurrence 
of equilibrium, periodic (including quasiperiodic) and 
chaotic regimes. Accordingly, the tools we employ are the 
time series and the spectra of Lyapunov exponents. The 
latter is particularly useful as a relatively sensitive tool 
to characterise the dynamics, with the Lyapunov spec- 
tra of the types (-, .), (0, -, -, . . .), (0, 0, -, . . .) 
and (+,0,—,...) corresponding to equilibrium, periodic, 
quasiperiodic (with two periods) and chaotic regimes re- 
spectively. Also to keep the numerical costs reasonable, 
the resolution of D in all the figures was, unless stated 
otherwise, taken to be D = 5. 



2 These authors seem to confine themselves to this value of 
v in order to obtain chaotic behaviour, for which one requires 
a > b + 1 (Sparrow 1982). This amounts to the expectation 
that a relaxes much more slowly than the magnetic field. 



Now given the fact that in many astrophysical settings 
(including that of the sun) the sign of the dynamo number 
is not known, we shall also study the effects of changes in 
the sign of D. 

We note also that the afl dynamo concept becomes 
invalid if D exceeds a certain limit (Choudhuri 1990). 
Furthermore, in general, as D is increased more modes 
(higher N) are required to achieve convergence (numeri- 
cally bounded solutions). 

3.1.1. Results for positive dynamo numbers 

For the sake of comparison with S&S91, we studied the 
dynamics of the system (15, 16, 23), for different values 
of the truncation order N. A summary of our numeri- 
cal results is given in Fig. 1 which is a plot of the two 
largest Lyapunov exponents as a function of the dynamo 
number for different truncations. In the following figures, 
the largest Lyapunov exponent is depicted by a solid line 
and its negative, zero and positive values indicate equi- 
librium, periodic and chaotic regimes. The simultaneous 
vanishing of the second Lyapunov exponent would imply 
the presence of quasiperiodic motion with two frequencies 
(i.e. motion on a 2-torus). It was not necessary to plot the 
third exponent, since no motion on T 3 or higher dimen- 
sional tori was observed which is not surprising in view of 
the results of Newhouse et al. (1978). 



Dynamo number D 



Fig. 1. Graphs of the two largest Lyapunov exponents for 
N = 2, 4, 6, 7, 8 and 10 (increasing downwards) for the case 
where / oc J ■ B and D > 



For a more transparent comparison, we have also pro- 
duced in Fig. 2 an analogous figure for the system con- 
sidered in S&S91. As the comparison of the Figs. 1 and 2 
shows, the main differences produced by the replacement 
of A^B^ by J • B are as follows: 

1. The chaotic regimes become less likely in the J ■ B 
case, in the sense that the intervals of the dynamo 
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Fig. 2. Graphs of the two largest Lyapunov exponents for 
N = 2, 4, 6, 7, and 8 (increasing downwards) for the case 
where / oc A^B^, and D > 




0.04 0.05 0.06 0.07 0.08 0.09 1 

Magnitude of the initial vector 



Fig. 3. Fragility in the dynamics with respect to small changes 
in the dynamo number D and the magnitude of the initial 
vector (A„, B„,C„). A cross represents a fixed point while a 
circle represents a limit cycle. The initial conditions correspond 
to randomly chosen vectors of specified magnitude 



number D over which the system is chaotic decrease 
dramatically. 

There exist indications for the presence of "multiple 
attractors" over substantial intervals of D, consisting 
of equilibrium and periodic states. These can be seen 
as regions of spiky behaviour in the solid line in Fig. 
1, for certain truncations (TV = 4, 6, 7, 8, 10). The be- 
haviour of the system alternates between fixed point 
solutions (where all exponents are negative) and pe- 
riodic orbits (where only the first one is zero) as the 
dynamo number D is slightly changed. 
The presence of such behaviour is potentially of great 
interest since it suggests that there exist intervals of 
D in which small changes in D can drastically change 
the behaviour of the system. This is also interesting, 
if one considers settings in which D or the initial con- 
ditions (IC) can vary slightly, but randomly, as the 
resulting behaviour would look very much like inter- 
mittency. To highlight this we have plotted in Fig. 3 
the behaviour of the TV = 4 truncation as a function of 
small changes in the dynamo number and the IC. As 
can be seen, small changes in cither D or IC can pro- 
duce important changes in the behaviour of the system. 
This therefore shows that there are substantial regions 
of D over which the behaviour of the system is sensi- 
tive to small changes in D and IC. Further, we have 
checked that this fragility is itself robust in the sense 
that taking a finer mesh of D does not qualitatively 
change this overall behaviour. 

Regarding the overall behaviour of the systems with 
respect to increases in TV, we observe the following. 
For small dynamo numbers, the behaviour seems to 
settle down to equilibrium and periodic states as TV is 
increased. For example as can be seen from Fig. 1, for 
dynamo numbers up to D w 900, the behaviour settles 
down for TV > 7. For larger values of D, however, we 



observe an increase in the dominance of the "multiple 
attractor" regime for the values of TV considered here. 
It is likely, however, that with increasing TV, these in- 
tervals only establish themselves at higher values of 
D. 

4. The transition to chaos appears to be very abrupt in 
the TV — 2 case, with the system going from a fixed 
point into a chaotic regime very rapidly, at least to 
within a resolution of AD w 10~ 4 , with no intermedi- 
ate behaviour being observed. For the case TV = 3 the 
system goes from a fixed point — > limit cycle — > chaos. 
For still higher TV, our calculations indicate that chaos 
becomes scarce. 

5. Chaotic regions were also found in the "multiple at- 
tractors" region, which were fragile with respect to 
small changes in the IC and the choice of D. 

3.1.2. Results for negative dynamo numbers 

Our results for the negative dynamo numbers are shown 
in Fig. 4. Also, in view of the sparseness of the results 
reported in S&S91 for the models with negative dynamo 
numbers, we present Fig. 5 as an analogous figure for their 
case. 

The main features of these models are: 

1. The chaotic regimes seem to become less likely in the 
A^B^ case. In fact, for the mesh size in D taken here, 
we only observed chaotic solutions in the case of TV = 8 
and then only for very high dynamo numbers. 

2. There are substantial intervals (in D) of "multiple 
attractors" (consisting of equilibrium and periodic 
states) for the A^B^ case. 

3. In both cases the behaviour for D <; —900 stabilises 
as TV is increased. This occurs for TV > 2 for the equi- 
librium regime and TV > 7 for periodic regime. These 
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Dynamo number D 

Fig. 4. Graphs of the two largest Lyapunov exponents for 
N = 7, 8, and 10 (increasing downwards) for the case when 
/ cx J • B and D < 




; ! 



Dynamo number D 

Fig. 5. Graphs of the two largest Lyapunov exponents for 
N = 4, 6, 7, and 8 (increasing downwards) for the case when 
/ cx A^B^ and D < 



results also indicate that there are parallels between 
the A^B^, case with negative dynamo numbers and 
the J ■ B case with positive dynamo numbers. In both 
cases, multiple attractor regions seem to dominate for 
large D values, as N is increased. 
4. For high N(> 7) in the J ■ B case, the transition is 
from T 2 to chaotic behaviour. This does not seem not 
true for N = 4,5,6 where the chaotic behaviour seems 
to appear abruptly. 

3.2. Case (II): f = Wi J ■ B + W 2 a\B\ 2 , with Wi ^ and 
W 2 ^0 

To study the effects of including the a|£?| 2 term, we 
use the dynamic a equation from Kleeorin & Ruzmaikin 
(1982) without a damping term 3 proportional to l/T 



The inclusion of this term will be studied in Sect. 4. 



da.M _ vt d 2 a M 
dt ~ ~B? dx 2 
1 



(29) 



Mop 



(V x £?)•£? — 



ao cos x — <xm 



B\ 



where /? is the combined (turbulent plus ohmic) diffusion 
of the field, p the density of the medium and po the mag- 
netic constant. 

Now using expression (2) for B and turning the system 
in a non-dimensional form using the same transformations 
as before, we obtain 



\B 2 = 



(30) 



where Ti — a is a dimcnsionless constant. This allows 
the analogue of the Eq. (11) to be written in the form 



dC d 2 C 



dt 



dx 2 



dA8B d 2 A B 
dx dx dx 2 



dA 



(31) 



+T 3 x {Dcosx-C) | B z +T 1 [ — 



where T 2 — ^ B ° and T3 = -. 



R A Bt, 



are dimcnsionless 



constants. Considering the same boundary conditions and 
spectral expansions as in the J • B case, Eq. (31) becomes 



C>Cn 

dt 



(N N \ 
^^H{n,m,l)A m Bi\ (32) 
m=l 1=1 J 

N N 

+r3^^DB ro B,Fi(n,m,l) 

m=l 1=1 
N N 

+r 3 ri ^mlDA m AiH 2 {n,m,l) 

m=l 1=1 

N N N 

C m BiB k H 3 (n, m, I, k) 

m=l 1=1 k=l 
N N N 

-T3T1 J2J2 lkCmAiAkHi ( n > m > l > k ^ 



m=l 1=1 k=l 

where i?'s are given by 
1 



Hi(n,m,l) = — 



m + 1 



7r \ (m + 1 + I — n) (m + 1 — I + n) 
m + 1 

(to + 1 + I + n) (m + 1 — I — n) 
to — 1 

(to — 1 + I — n)(m — 1 — I + n) 



(33) 
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TO 



H 2 (n,m,l) = - 



(to — 1 + I + n) (to — 1 — I — n) J 
1 ( n-l 



7r \ (n — I + m + l)(n — I — m — 1) 
n — I 



(34) 



(n — I + to — l)(n — I — to + 1) 
n + Z 

(n + I + to + l)(n + Z - to - 1) 
n + Z 

(n + l + m- l)(n + Z - to+ 1) 



H 3 (n, to, Z, fc) 
H±{n, to, Z, fc) 



+ 
+ 
+ 

- [<5(to — Z, fc — n) — <5(m — l,k + n) (35) 
— <5(to + Z, fc — n) + (5(m + Z, fc + n)] , 

i [5(m -n,l-k) + S(m -n,l + k) (36) 
—S(m + n,l — k) — S(m + n, I + k)] . 



Note that S(n, to) is 1 if n — m = but 2 if n = m = 
and Hi = H 2 — if n + m + Z + 1 is even. 



3.2.1. Results 

Our results of the study of the system (15), (16) and (32) 
for positive dynamo numbers are depicted in Table 1. As 
can be seen, the effect of the inclusion of the a|-B| 2 term is 
dramatic and seems to eliminate the possibility of chaotic 
behaviour for all N. 



Table 1. Results for the case (II) for D > 0. D\ indicates 
where the origin becomes unstable as a fixed point and D2 the 
dynamo number where all fixed points become unstable and 
the solution becomes periodic 



N 


Di 


D 2 


2 


10 


> 2000 


3 


15 


> 2000 


4 


115 


115 


5 


205 


205 


6 


205 


240 


7 


235 


240 


8 


250 


250 



For the lower truncations of N = 2 and 3, we only 
observe fixed point solutions for all D up to D w 2000. 
For higher order truncations, with moderate D, there is a 
sequence of fixed points followed by stable periodic cycles. 

The corresponding results for the negative dynamo 
numbers are shown in the Table 2, and again this is very 
similar to Table 1 with no evidence for chaotic behaviour 
at small and moderate D. In this case the N = 2 system 
has the origin as the fixed point for D down to —2000. 



Table 2. Results for the case (II) for D < 0. Di indicates 
where the origin becomes unstable as a fixed point and D 2 the 
dynamo number where all fixed points become unstable and 
the solution becomes periodic 



N 



D 2 



2 


< -2000 


< -2000 


3 


-70 


-80 


4 


-85 


-95 


5 


-85 


-95 


6 


« -95 


-100 


7 


« -95 


^ -100 


8 


w -95 


^ -100 



4. Case (III): Robustness with respect to changes 
in the damping term 

In this section we employ the equation proposed by Klee- 
orin et al. (1995) in the form 



(37) 



da M a M 
dt ~ T 

-^( { y xB )-B- aaCOSX - aM \B\ 2 
Mop V P 

as the evolutionary equation for the back reaction of the 
magnetic field on the time dependent part of a. In the 
above equation T is the characteristic time on which the 
small scale magnetic helicity changes, which is typically 
much longer than the turbulent diffusion time scale. 

Using the same expression for B from Eq. (2) and 
proceeding in the same way as in the previous cases we 
obtain the differential equations for C n to be 



dCn 

dt 



N N 



-TiC n + r 2 E H ( n > m > ^ AmBl 

\m=l 1=1 / 

N N 

m=l 1 = 1 
N N 

+T 3 Ti ^2mlDA m AiH 2 (n,m,l) 

m=l 1=1 
N N N 

-r 3 E E C m BiB k H 3 {n, to, Z, k) 

m=l 1=1 k=l 
N N N 

-r 3 ri Y ^2^2lkC m AiA k H A (n,m,l,k), 



(38) 



m=l 1=1 k=l 

where T 4 = is a dimcnsionless constant. 
4-1. Results 

Our results of the study of the system (15), (16) and (38) 
are shown in Tables 3 and 4. Although more modes are 
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required in order to obtain convergence for higher dynamo 
numbers, the results shown in Table 3 and 4 seem to in- 
dicate that this type of change in the damping term does 
not produce qualitative changes in the behaviour of the 
system. This is reasonable, since the functional forms of 
the modal equations are quite similar in Eqs. (32) and 
(38). 

The inclusion of the cxm/T term does not change the 
qualitative behaviour of the smaller truncations (N = 2 
and 3 for D > and N = 2 for D < 0, where we observe 
only fixed points as before). At moderate dynamo num- 
bers, the qualitative behaviour is almost the same and re- 
mains periodic for \D\ > \D 2 \, but D 2 is changed slightly. 



Table 3. Results for the case (III) for D > 0. Di indicates 
where the origin becomes unstable as a fixed point and D 2 the 
dynamo number where all fixed points become unstable and 
the solution is a periodic orbit 



N 


Di 


D 2 


2 


10 


> 2000 


3 


15 


> 2000 


4 


110 


115 


5 


170 


170 


6 


« 200 


?a 200 


7 


w 200 


« 200 


8 


w 200 


« 200 



Table 4. Results for the case (III) for D < 0. D\ indicates 
where the origin becomes unstable as a fixed point and D2 the 
dynamo number where all fixed points become unstable and 
the solution is a periodic orbit 



N 


Di 


D 2 


2 


< 2000 


< 2000 


3 


-70 


-80 


4 


-80 


-95 


5 


-95 


-95 


6 


« -95 


« -105 


7 


« -95 


« -110 


8 


« -95 


« -110 



We also note that all systems considered here, in par- 
ticular Cases (II) and (III), have a common pattern of 
behaviour, namely that as D is increased, A and B oscil- 
late with slowly increasing amplitudes about zero. On the 
other hand, qm oscillates with an increasing amplitude 
around a rapidly increasing average. Also if D > 0, olm 
oscillates about a positive average and about a negative 
average for D < 0. 



5. Conclusions 

We have studied the robustness of truncated afl dynamos 
including a dynamic a equation, with respect to physically 
motivated changes in the driving term and a change in 
the damping term appearing in the dynamical a equation. 
We studied these systems with respect to changes in the 
dynamo number D, the truncation order N and the IC. 
Our results show that the changes in the driving term 
have important effects on the dynamical behaviour of the 
resulting systems. In particular we find that 

— chaos is much less likely in systems with a driving term 
of the form J B (with positive D), as opposed to those 
involving A^B^. 

— the inclusion of the a|-B| 2 term has a dramatic effect in 
that it suppresses the possibility of chaotic behaviour 
at moderate dynamo numbers. 

— changes in the sign of the dynamo number can also pro- 
duce important changes. In the case where the driving 
term is given by A^B^, using D < makes chaotic 
behaviour much less likely (which seems to be the mir- 
ror image of the case where the driving term given by 
J ■ B and D > 0). 

— in case (I) there exists substantial intervals of D for 
which the systems seem to possess "multiple attrac- 
tors" (consisting of equilibrium and periodic states). 
As a result small changes in either D or the IC can 
produce important changes in these regimes. This form 
of fragility can be of importance, especially in presence 
of noise, where the system would behave in an inter- 
mittent way. 

Finally to recapitulate our motivation for studying dif- 
ferent formulations of dynamic a feedback, we note that 
even the usual expression for the driving term, J B, de- 
rived from first principles could still be inappropriate, as it 
involves uncontrolled approximations. However, it is clear 
that / has to be a pseudo-scalar (because a is a pseudo- 
scalar), and the most obvious possibilities are indeed the 
ones that we have studied. Our investigations have shown 
that the actual choice can significantly alter the over- 
all conclusion. Therefore, all conclusions, especially those 
concerning the occurrence of chaos, should be taken with 
utmost care. 
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